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Abstract 

The thermomagnetic convection of magnetic fluids in a cylindrical geometry subjected 
to a homogeneous magnetic field is studied. The study is motivated by a novel thermal 
instability [W. Luo et al, Phys. Rev. Lett. 82, 4134 (1999)]. As model system a composite 
cylinder with inner heating is considered which reflects the symmetry of the experimentally 
setup. The general condition for the existence of a potentially unstable stratification in the 
magnetic fluid is derived. Within a linear stability analysis the critical external induction 
for the onset of thermomagnetic convection is determined for dilute and nondilute magnetic 
fluids. The difference between both thresholds allows to test experimentally whether a test 
sample is a dilute fluid or not. 
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1 Introduction 



Magnetic fluids (MFs) are superparamagnetic fluids formed by a stable colloidal suspension of 
ferrimagnetic nanoparticles dispersed in a carrier liquid The behaviour of MFs is character- 
ized by the complex interaction of their hydrodynamic and magnetic properties with external 
forces. For the phenomenon of thermomagnetic convection these forces are a temperature gra- 
dient and a magnetic field. Most studies consider the geometry of a horizontal layer which is 
simultaneously subjected to a vertical temperature gradient and either to a constant vertical 
magnetic field ||, |3], | or to a vertical magnetic field with a constant gradient ||. Studies 
for a cylindrical geometry are few || |7| and focus on the thermomagnetic convection under 
microgravity ||. 

The present analysis of thermomagnetic convection in a cylindrical geometry is motivated 
by a recently observed novel convective instability. In ||, 10] a horizontal layer of MF between 



two glass plates is locally heated by a focused laser beam. It passes perpendicularly through the 
layer in the presence of a homogeneous vertical magnetic field. The absorption of the light by 
the fluid generates a temperature gradient and subsequently a refractive index gradient. This 
gradient is optically equivalent to a diverging lens, leading to an enhancement of the beam 
divergence. As result, a stationary diffraction pattern of concentric rings is observed for zero 
magnetic field. Above a certain threshold of the magnetic field, the circular rings are replaced 
by polygonally shaped patterns which switch among different shapes alternatively. Based on 
numerics it was stated that the characteristic time scales for mass and thermal diffusion are 
equal || . Thus thermal conduction and thermal diffusion contribute equally to the formation of 
the diverging lens. The polygonally shaped diffraction patterns were interpreted as 'fingerprints' 
of vertical convection columns [|l0|] . 

Both claims are controversial due to the following reasons. An analysis of the physical 
quantities reveals that both characteristic time scales are several orders of magnitude apart 



(below and |jT| ). Direct spatial temperature measurements [12] and independent measurements 



of all relevant material parameters for organic dispersions showed that the ring pattern is 



essentially caused by the temperature contribution and not by the concentration contribution. 
The determination of the relevant time scales confirmed that the experiment is dominated by 
the characteristic time for convection. Presently there is no sound theoretical description for 
the dependence of the size of the outmost of the concentric rings on the laser power (see Fig. 2 
in m 



A major hindrance of the studied system in ||, 1C] is that it is almost impossible to gain 
information about the spatial distribution of temperature and concentration inside the MF 
layer. Due to the lack of sound internal information many hypotheses can be brought forward 
to explain the experimental results. Among them is the recent discussion whether the shape 
instability of a hot nonmagnetic bubble surrounded by MF can be accounted for the observed 



phenomena [14, 15 1. 

This situation motivates the present work in which a model system is studied which reflects 
the essentials of the experimental setup. These are the axis-symmetry of the heating and the 
finite hight and width of the layer. The aim is to determine the necessary conditions and the 
critical magnetic inductions for the appearance of axial convection columns. Thus the focus 
is on the general requirements for thermomagnetic convection in an axis-symmetric heated 
system. 

This paper is organized as follows: The system and the relevant equations of the problem 
as well as the condition for a potentially unstable stratification in the fluid are displayed in 
the next section. Based on a linear stability analysis (Sec. Ill), the results are presented and 
discussed in Sec. IV. In the final section, the results are summarized. 



2 Model and Equations 

The model system is given by a composite, circular cylinder of height h which consists of three 
parts. The inner cylinder of radius Ri, constant temperature T%, and constant susceptibility 
x(?l) is surrounded by a middle cylinder of radius i?2- In the gap i?2 — Ri the temperature 
decreases to To < T% and consequently \ is a spatially varying quantity, % = x[T(r)]. The 
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outer cylinder has the radius i? ut j where in the region i? ut — R2 the constant temperature To 
and constant susceptibility x(Tq) is present. The whole system is subjected to a homogeneous 
vertical magnetic field and its effective susceptibility is given by 

-R 2 



Xcs = -^—(rix(Ti) + (R2 - Ri) [ 2 dr X [T(r)] + (R out - R 2 )x(T )) . 

-n-out V JR 1 J 



(1) 



In the presence of a uniform external magnetic induction B ex t, the internal field in the gap 
is given by H int = B ext /[/x (l + Nx)]. The susceptibility of the MF is x = XlO- + PiXl), 
where xl is the susceptibility according to Langevins theory which assumes non-interacting 
particles. Higher order terms in xl are included in order to determine the magnetic (or Kelvin) 
force density beyond the dilute limit X = XL- The coefficient f3\ was determined in different 
microscopic models ]16| , |17| , [18[ ] which all provide the same value (5\ = 1/3. The demagnetization 
factor N accounts for the finite size of the composite cylinder and is a function of the heigh-to- 



diameter ratio 7 = h/(2R ont ) and the effective susceptibility Xeff llMl- The Kelvin force follows 



then as J20] 



* K = — ' W 

XL 



where 



Xl {N + ft [3N XL (1 + 01XL) ~ 1]} 
(l + iV X )3 



jp Ai l J " ~ Ml I'-' 1 y AL \-l ~ yLALJ J-J/ /Q x 
b Xh = , AT- 7\3 • i 6 ) 



Considering MFs as binary mixtures, it is necessary to evaluate the influence of temperature 
and concentration on pattern phenomena by analyzing the corresponding time scales. These 
are the characteristic time for convection t c = L\j k and for mass diffusion = L^/D, where 
L c (Ld) is the typical length for convection (diffusion), k the thermal diffusivity, and D the 
mass diffusion coefficient. Using the data given in [Q, [|, one gets k ~ 4 x 10 -8 m 2 s _1 and 
D ~ 8 x 10~ 12 m 2 s _1 . Except in special designed geometries as in [21|, where L c ^> L^, the 



characteristic time for diffusion is of three orders of magnitude larger than the characteristic 
time for convection. Since in our model both length scales are equal to the gap width R2 — Ri, 
diffusion phenomena can be neglected. 



The system is governed by the equation of continuity, the Navier-Stokes equations, and the 
equation of heat conduction for the MF which are in nondimensional form 

divv = 0, (4) 

-_ + (vgrad)v = V (-gradp + Av) + MF XL ^-^ — , (5) 
df 

_ + (vgrad)T = AT, (6) 

where the Prandtl number V = v / 'k characterizes the fluid and the magnetization number 
M = i?c xt (i?2 — Ri) 2 /{l^opn 2 ) tunes the external excitation. Denoting v as kinematic viscosity, 
the velocity v = (u, v ) is scaled with ft/(i?2 — Ri), time with (R2 — R\) 2 /k, temperature with 
(Ti — To), and pressure p with pnv /(R2 — Ri) 2 ■ A and grad are the corresponding differential 
operators in the plane cylindrical coordinates r and 4>. Rigid boundary conditions are assumed 
for the velocity at the inner and outer radius of the gap, u = dfU = at f = — if) and 
f = 1/(1 — 77), where the radii ratio is given by 77 = i?i/i?2- The temperature is assumed to be 
constant at each boundary, T(f = 77/(1 — 77)) = T\ and T(f = 1/(1 — 77)) = To. 

Since the Kelvin force is the only destabilizing force present in the system, one has to 
determine which profile leads to a potentially unstable stratification in the fluid. For heating 
at the inner radius, the required profile is given in Fig. [l|(a): the r-component of the Kelvin 
force density has to act inwards and its absolute value increases monotonically outward. With 
such a profile a fluid volume at the distance r + 5r (solid rectangle in Fig. |l](b)) experiences a 
larger force towards the center compared to a fluid volume at the distance r (dashed rectangle) . 
Moving the latter fluid volume from r to r+5r (dot-dashed rectangle) results in an effective force 
which points in the direction of the displacement (indicated symbolically by the subtraction of 
the arrows in Fig. 0(b) bottom). This force may enhance small displacements of warmer fluid 
volumes towards cooler regions and thus making the stratification potentially unstable. 

The above argument has to be tested for the quiescent conductive state which is given by 
vg* = and Tq = Tq + (T\ — Tq) ln[f(l — 77)]/ In 77. Applying the condition for a destabilizing 
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force profile to the r-component of the Kelvin force density in Eq. (g) leads to the condition 

for all r £ [77/(1 — 77), 1/(1 — 77)] . (7) 



Since df^K,r is a monotonously decreasing function of f, it is sufficient if df^K,f < at 
f = rj/(l —77) in order to fulfil the condition ([?]). Depending on the temperatures T\ and 
To, the Langevin susceptibility xl, the demagnetization factor N, and j3\ the condition (|7|) 
entails that the radii ratio 77 has to be larger than a critical value. For realistic temperatures 
T\ above a room temperature of To = 300 K, it becomes clear that this condition is met 
only in a narrow gap (see Fig. |2|). This is plausible because the nonlinear temperature profile 
Tq ~ ln[f(l — 77)] can be well approximated in a small gap by a linear profile which always 
satisfies the requirement (Q) if /3i = 0. 

3 Linear Stability Analysis 

Exploiting the smallness of the gap, in the linear stability analysis terms as df{df + 1/f) 
are approximated by d 2 and the new variable Q = f — 77/(1 — 77) is introduced. All small 
disturbances from the ground state are decomposed into normal modes, i.e. into components 
of the form [u,p,T] = e nt cos(l(p)[u(C),p(C),T(()] and v = e nt sm(l(f))v(C), respectively. The 
nondimensional growth rate is denoted by n and I is the azimuthal wave number. For marginal 
stability, n = 0, the differential equations to solve are 

d 2 2 V_ a 2 ( d 2 2 \_ 2 M £ f df G 



dc 2 a r--p\w~ a r = ~ a ^ hL w^' (8) 



where a = (1 — rj)l/rj and 

+4N X 2 MN -4^) + Xl(N - 10/30 + 2 A - 2N] . (10) 
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In order to satisfy the four boundary conditions u = d^u = at £ = 0, 1, the ansatz 

K 

u(C) = X] Pm sinh(aC) + B m cosh(aC) + C m sin(&C) + D m cos(&C) (11) 
m=l 

with a = y / q 2 n + a 2 and 6 = \f q 2 ^ — a 2 is chosen. g m is a root of a transcendental equation 
and the constants B m , C m , and D m are determined by the boundary conditions (for details see 
H). With Eq. © the solution of Eq. (|) reads f(() = £* =1 P m [T m {Q + Cie°< + C 2 e"^] . 
T m (£) is the solution of the inhomogeneous equation (due to its lengths not given here) and 
the constants C\ and C2 are determined by the boundary conditions T = at £ = 0, 1. Using 
Eq. (|ll|) and the solution for T(£), Eq. @ can be approximately solved by the Galerkin method. 
Due to the good convergence, all presented results are based on the third approximation (see 
Table p]). For the calculations fluid parameters of EMG 901 are used: p = 1.53 x 10 3 kg m , 
v = 6.54 x 10~ 6 m 2 s _1 , xl = 3, and k = 4.2 x 10" 8 m 2 s _1 The temperature at the 

outer radius of the gap R2 = 1 cm is fixed at To = 300 K. The hight of the composite cylinder 
is given by h = 1 cm and the inner radius by rj = 1.01r/ c (/?i = 1/3, N = 1). The choice of 
VciPi = 1/3, N = 1) (solid line in Fig. ^) ensures that for all following parameter sets the 
condition (^) is fulfilled. 



4 Results and Discussion 

Solving Eq. @ with the Galerkin method and subsequently minimization with respect to the 
azimuthal wave number determines the critical external induction B c and the corresponding 
wave number l c (Figs. || and ||). Four different parameter sets were chosen: a dilute {(3\ = 0) 
and a nondilute {(3\ = 1/3) MF with i? ut = 00 {N = 1) and -Rout — 3.33 cm, respectively. The 
demagnetization factor N for the resulting height-to-diameter ratio 7 = 0.15 and the effective 
susceptibility XcS of the composite cylinder accordingly to Eq. (||) is taken from the data given 
in®. 

Decreasing the temperature difference from AT = 70 K to AT = 4 K causes a dramatic 
increase in the critical induction of nearly three orders of magnitude (Fig. |3|). With decreasing 
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temperature difference the critical radii ratio grows, i.e. the allowed gap becomes more narrow. 
Since the convection rolls prefer the same length scale in r and (^-direction, much more rolls 
have to be driven in a very small gap. The energy for this effort comes from the external 
induction which is why it amplifies drastically for small AT (Fig. 

Whereas the critical azimuthal wave number l c is independent of x{xl) and N (Fig. ||), the 
critical induction varies. First the two thresholds for the case of an infinitely extended layer, 
-Rout = oo, are compared. The inclusion of a quadratic term in the susceptibility with (3\ = 1/3 
results in a lower threshold for the onset of convection than in the dilute case (3\ = (solid and 
long-dashed line in Fig. ||). The difference between the thresholds is nearly the same value, 
B c (f3 1= l/3,N=1) ~ 0.63£ c (/3i=0,iV=l), for all tested temperatures 304 K < T x < 370 K. 

Now the thresholds for the case of a finite layer with R ou t — 3.33 cm are compared. Contrary 
to the previous case, the threshold for a dilute fluid (dot-dashed line in Fig. ||) is lower than 
for a nondilute fluid (dotted line). Again the difference is almost constant over the entire 
temperature range, 5 c [/?i=l/3, iV( 7 =0.15, Xeff)] - 1.18S c [/3i=0, iV( 7 =0.15, Xeff)]- 

The relation of the different thresholds is caused by the value of f XL for the given com- 
binations of N, Pi, and \l- f\L can be considered as a measure for the strength of the 
magnetic force in the gap: as higher the value of f XL as lower the critical external induction 
necessary to trigger the convection. Figure [5] shows the value of f XL for the four considered 
parameter sets. At \L = 3 the relation f XL ((3i = 0, N = 1) < / XL (/?i = l/3, N = l) (see cross- 
sections of the long-dashed and the solid line with the vertical solid line) is the reason that 
the threshold for the dilute fluid is higher than for the nondilute fluid. The opposite relation 
f XL [Pi=0, iV(7=0.15, Xefr)] > f\L [/^i=l/3, AT (7=0.15, Xeff )] (see cross-sections of the dot-dashed 
and the dotted line with the vertical solid line) causes the opposite relation for the thresholds 
in the case of a finite layer. 

The physical reasons which cause theses differences are the following. The Kelvin force is 
proportional to the magnetization in the magnetic fluid. Thus as higher the magnetization 
is, as lower the external induction can be in order to generate the same strength of the mag- 



S 



netic force. In the infinite case, where N = 1 is independent of x> a higher concentration of 
magnetic particles in the fluid leads to a higher magnetization and therefore to a lower thresh- 



old. In the finite case, the demagnetization factor depends on x [19] and is smaller for higher 
concentrations of magnetic particles than for lower concentrations. A higher demagnetization 
factor means a lower inner field and a higher magnetization, respectively. Therefore in the 
finite case an increase in the concentration results in two counteracting effects with respect 
to the magnetization. In the studied example of xl = 3 the influence of the demagnetization 
effect wins: the dilute fluid has the lower threshold. But for xl > 3.2 the direct influence of 
the concentration succeeds over the demagnetization effect. The nondilute fluid has the lower 
threshold (the dotted line is then above the dot-dashed line, see Fig. ||). 

In the infinite and finite case the clear and measurable difference between the thresholds 
opens a very good opportunity to decide whether a test sample is a dilute fluid or not. Just 
by measuring the threshold for the onset of convection in the proposed model system the 
answer can be given. The critical induction depends on the fluid and system parameters as 
B c ~ Ky / p/(i?2 — R\)- By choosing fluids with low (high) density and thermal conductivity 
and a large (small) radius R2, the threshold can be lowered (raised) corresponding to the 
experimentally available magnetic fields. 

The major obstacle in order to compare the results with the experimental data in JlOj j is 
the lack of an experimentally determined spatial profile of the temperature inside the sample. 
Therefore it is not possible to extract an estimation what might be the values of Ri and R2 
in the experiment. Nevertheless the calculated values indicate that really high critical exter- 
nal inductions are necessary to trigger vertical convections rolls by a pure radial temperature 
gradient. The threshold for the induction reaches extremely high values if one extrapolates 
towards radii in the range of hundreds of micrometers, not unlikely due to the focused laser 
beam used in the experiment |], |l(]] . This leads to the conclusion that vertical convections rolls 
due to a pure radial temperature gradient are unlikely to account for the observed phenomena 
of polygonally shaped diffraction patterns. 
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Due to the lack of information from inside the sample, it is not clear whether the temperature 



profile in the experimental sample is purely radial. There are hints in 24] that due to the 
absorption along the way of the laser beam a vertical temperature distribution exists as well. 
A further cause for an axial temperature gradient is the heat loss through the glass plates by 
which the MF layer is sandwiched. If such a vertical temperature gradient comes into play, 
concentration gradients due to the Soret effect may become important. The relative influence of 
temperature and concentration gradients is strongly effected by the relation of the characteristic 
times t c and which depend quadratically on the lengths L c and L^, respectively. With a 
radial and a vertical temperature gradient present, it becomes even more important to have 
reliable data of the internal profiles to estimate these lengths. 

Two differences between the model and the motivating experiment should be noted. The 
constant temperature T% for the inner cylinder is not given in the experiment. How much the 
temperature varies in this inner area is not known. The numerical calculations in || suggest 
a difference of about 15 K. The thresholds were calculated for rigid boundaries whereas in 
the experiment the fluid layer boundary is free. For thermal convection in a rotating layer of 
magnetic fluid the influence of rigid and free boundaries on the threshold was studied in Q. 
Considering the case of zero rotation, the thresholds differ by not more than 20% (see Figs. 3 
and 4 in M for rotation number T — > 0). With respect to the above mentioned major obstacle, 
these differences may alter the results only marginally. 



5 Summary 

A model system of a composite cylinder of finite size with axis-symmetrical temperature dis- 
tribution is presented. Using the Kelvin force density (|2],[|) a general condition (0) is derived 
for which a potentially unstable stratification exists if the inner cylinder is heated. Depending 
on the temperature difference, the size of the composite cylinder and the dilute or nondilute 
character of the magnetic fluid, the critical gap sizes are calculated. The general result is that 
only in a narrow gap the requirement for a potentially unstable stratification is met. Exploiting 
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this property, a linear stability analysis is performed in order to determine the critical external 
induction for the onset of thermomagnetic convection. With decreasing temperature difference, 
the critical induction increases dramatically. The reason is that for smaller temperature differ- 
ences which demand smaller gaps in order to fulfil condition (0), much more convection rolls 
have to be driven. The driving of these many rolls causes the drastic increasing of the critical 
induction. The distinct difference between the threshold for dilute and nondilute magnetic 
fluids allows to use the considered system for an experimental determination whether a test 
fluid is a dilute or nondilute one. 

The rather high external induction, needed to stimulate the convection flow, leads to the 
conclusion that vertical convections rolls due to a pure radial temperature gradient are unlikely 
to account for the observed diffraction patterns. The consideration of a vertical temperature 
gradient entails that concentration gradients may become relevant. To answer this question the 
characteristic time for the diffusion with respect to whose for convection has to be re-estimated. 
Because they will not be necessarily apart by orders of magnitude as in the case of a pure radial 
temperature gradient. Also different convection patterns can be expected with the presence of 
a vertical temperature gradient. In order to come to a correct statement about the contribution 
of mass and thermal diffusion to the diffraction patterns in magnetic fluids, spatial temperature 



measurements and independent measurements of all relevant material parameters as in [12, 13] 
for an organic dispersion are highly desirable. 
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Table 1: Critical external induction B c in dependence of the order of approximation for N = 1. 
T^K] ft B c [T] 

K=l K=2 K=3 K=4 
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0.02663 
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Figure Captions 

Figure 1: Required destabilizing force profile of the radial component of the magnetic force 
density fx, r for inner heating (a). A fluid volume at the distance r + 5r (solid rectangle, (b)) 
experiences a larger force than a fluid volume at the distance r (dashed rectangle, (b)). Moving 
the latter volume from r to r + 5r (dot-dashed rectangle, (b)) results in an effective force which 
points in the direction of the displacement (dot-dashed arrow, (b)). 

Figure 2: Region of potentially unstable and stable force profiles for a fixed outer temper- 
ature of To = 300 K and xl = 3. The four different sets are ft=0, N=l (long-dashed line), 
/3i=0, iV=0.7 (dot-dashed line), ft=l/3, N=l (solid line), and ft=l/3, iV=0.7 (dotted line), 
where the first and last one practically coincide. 

Figure 3: Critical external induction B c versus inner temperature Tl for a room temperature 
of To = 300 K. For a horizontally infinitely extended layer, i.e. Rout = oo, the inclusion of a 
quadratic term in the susceptibility with ft =1/3 (solid line) results in a lower threshold for the 
onset of convection than in the dilute case, ft = (long-dashed line). Contrary for R out ~ 3.33 
cm, the critical induction for the dilute fluid (dot-dashed line) is lower than for the nondilute 
fluid (dotted line). The fluid parameters of the magnetic fluid EMG 901 and the size of the 
cylinder are given in the text. 

Figure 4: Critical azimuthal wave number l c versus inner temperature Tl for To = 300 K. 
With decreasing temperature difference the wave number, i.e. the number of convection rolls, 
increases dramatically from Z c = 14 for AT = 70 K to Z c = 315 for AT = 6 (•). 

Figure 5: f XL , a measure for the strength of the magnetic force, versus the Langevin 
susceptibility xl- The four different sets are ft=0, N=l (long-dashed line), ft = 1/3, N=l 
(solid line), ft=0, N(j=0.15, Xcs) (dot-dashed line), and ft = 1/3, iV(7=0.15, Xcs) (dotted 
line). The vertical solid line at xl = 3 is a guide for the eye. 
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Figure 1: Lange, Physics of Fluids 




Figure 2: Lange, Physics of Fluids 
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Figure 5: Lange, Physics of Fluids 
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